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Abstract 

CIh' 

This study is aimed at making a calculation about the impact of the two most commonly used solar EUV flux 
. models - SOLAR2000 (S2K) of Tobiska (2004) and EUVAC model of Richards et al. (1994) - on photoelectron fluxes, 

1^, volume emission rates, ion densities and CO Cameron and CO2 UV doublet band dayglow emissions on Mars in 

O ' three solar activity conditions: minimum, moderate, and maximum. Calculated limb intensities profiles are compared 

' with SPICAM/Mars Express and Mariner observations. Analytical yield spectrum (AYS) approach has been used to 

. calculate photoelectron fluxes in Martian upper atmosphere. Densities of prominent ions and CO molecule in excited 

triplet a'^n state are calculated using major ion-neutral reactions. Volume emission rates of CO Cameron and COj UV 
doublet bands have been calculated for different observations (Viking condition. Mariner and Mars Express SPICAM 
observations) on Mars. For the low solar activity condition, dayglow intensities calculated using the S2K model are 
, ~40% higher than those calculated using the EUVAC model. During high solar activity, due to the higher EUV fluxes 

' at wavelengths below 250 A in the EUVAC model, intensities calculated using EUVAC model are slightly higher (~20%) 

than those calculated using S2K model. Irrespective of the solar activity condition, production of Cameron band due 
to photodissociative excitation of CO2 is around 50% higher when S2K model is used. Altitude of peak limb brightness 
of CO Cameron and COj UV doublet band is found to be independent of solar EUV flux models. Calculated limb 
intensities of CO Cameron and CO^ UV doublet bands are on an average a factor of ~2 and ~1.5, respectively, higher 
than the SPICAM Mars Express observation, while they are consistent with the Mariner observations. 

• 1-^ - 

d ■ 1 Introduction state coj (x^n) gives emission in ultraviolet wavelengths 

- - ■ at 288.3 and 288.6 nm. Apart from these emissions, Fox- 

First observation of CO Cameron and C0+ UV dou- Duffenback-Berger band of C0+ (A^n^ ~ X^nj, fourth 
blet emissions in the Martian dayglow were made by the p^^j^j^^ ^^^^ g^.^^ negative band of C0+ (B - 

Mariner 6 and 7 flybys in 1969-1970 (Barth et al., 1971; x), and several atomic line emissions of carbon and oxy- 
Stewart, 1972). These observations provided an opportu- ^^^^^ ^^^^ ^1^^ recorded by Mariner 6, 7, and 9 

nity to study the Martian upper atmosphere in a greater (g^^^j^ ^g^^. gtewart, 1972; Stewart et al., 1972). 

detail. The CO Cameron band (a^B ~ X^Y.^- 180 - ^-^^^ j^^lp theoretical calculations and laboratory 
260 nm) system arises due to the transition from the measurements, Barth et al. (1971) proposed possible mech- 
excited triplet a^H state, which is the lowest of triplet ^^^^^^ dayglow emission observed in the Martian 

states, to the ground state X1S+ of CO. Doublet transition atmosphere. Maximum intensity of CO Cameron band and 
(S2S+ _ X^Vl) from excited G0\ (B^E,,) to the ground 
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UV doublet observed by Mariner 6 and 7 were 600 kR at 
~131 km and 35 kR at 148 km, respectively. 

Emissions from Cameron band and COj UV doublet 
bands were also observed in 1971-1972 by Mariner 9, the 
first spacecraft to orbit Mars. Stewart et al. (1972) ob- 
served a reduction in the intensity of Cameron band by a 
factor of 2.5 compared to Mariner 6 and 7 observations. 
They attributed this difference to the reduction in the so- 
lar activity and better calibration of Mariner 9 instrument. 
The observed maximum slant intensities of CO Cameron 
band were between 200 and 300 kR and averaged topside 
scale height for the same band was 17.5 km. Stewart (1972) 
also observed a good correlation between CO Cameron 
band intensity and solar F10.7 flux, which suggest that 
these emissions are controlled by the incident solar photon 
flux. 

Since the Mariner 6, 7, and 9 UV observations, SPI- 
CAM (spectroscopy for the Investigation of the Character- 
istics of the Atmosphere of Mars) on-board Mars Express 
(MEX) is the first instrument dedicated for the aeronom- 
ical studies of Mars. SPICAM has broaden our under- 
standing about the Martian dayglow. Emissions observed 
by SPICAM in UV dayglow are H Lyman-a emission at 
121.6 nm, the atomic oxygen emissions at 130.4 and 297.2 
nm, the Cameron band {a^H — AT^E^) and fourth posi- 
tive band (A^H - X^T.+) of CO, and ultraviolet doublet 
band (i?^I]+ — A^II) emissions of COj (cf. Leblanc et al., 
2006; Chaufray et al., 2008). These emission features are 
similar to those observed by Mariner 6, 7, and 9 but 
with better sensitivity, and spatial and temporal cover- 
age. SPICAM has observed these dayglow emissions on 
Mars throughout the Martian year and showed the effect 
of solar zenith angle (SZA), seasonal variation, and Mar- 
tian dust storms on the dayglow emissions (Leblanc et al.. 
2006, 2007; Shcmatovich et al. , 2008; Simon et al., 2009; 
Cox et al., 2010). SPICAM also provided first observation 
of N2 UV emissions in Martian dayglow (Leblanc et al., 
2006, 2007). N2 Vegard-Kaplan VK (0, 5) and (0, 6) band 
emissions at 260.4 nm and 276.2 nm, respectively, have 
been observed; N2 VK (0, 7) emission at 293.7 nm has also 
been reported, but it has large uncertainty (Leblanc et al., 
2007). The detailed model of N2 dayglow emissions on 
Mars is presented elsewhere (Jain and Bhardwaj, 2011). 

Several theoretical studies have been made for 
the dayglow emissions on Mars (Fox and Dalgarno, 



1979; Mantas and Hanson, 1979; Conway, 1981; 
Shematovich et al., 2008; Simon et al., 2009). First de- 
tailed calculation of dayglow emission on Mars was carried 
out by Fox and Dalgarno (1979). Calculated overhead in- 
tensities of CO Cameron and CO J UV doublet bands were 
49 kR and 12 kR, respectively, for the low solar activity 
condition similar to Viking landing (Fox and Dalgarno, 
1979). Simon et al. (2009) used Trans-Mars model to cal- 
culate limb intensities of Cameron and COj UV doublet 
emissions for low solar activity condition (similar to Viking 
landing) and compared them with SPICAM-observations. 
Their calculated intensities are higher by ~25% than the 
observation. Simon et al. (2009) had to reduce the Viking 
CO2 density in the model atmosphere by a factor of 3 to 
bring the altitude of peak emission in agreement with the 
observation. 

Seasonal effects on intensities of Cameron and UV dou- 
blet bands have been observed by SPICAM (Simon et al., 
2009; Coxet al., 2010). Cox et al. (2010) have presented 
a statistical analysis of Cameron band and UV doublet 
emissions, peak altitude of emissions, and ratios between 
UV doublet and Cameron band. Averaged peak emis- 
sion brightness (altitude of peak emission) observed by 
Cox et al. (2010) for CO Cameron and UV doublet bands 
are 118 ± 33 kR (121 ± 6.5 km) and 21.6 ± 7.2 kR (119.1 
± 7.0 km), respectively. They also presented observations 
for one particular season, solar longitude (Ls) = 90 to 180°, 
and compared observational data with model calculations 
based on Monte Carlo code, which has been used also by 
Shematovich et al. (2008) for the Martian dayglow studies. 

Modelling of CO Cameron and CO^ UV doublet day- 
glow emissions requires a sophisticated input solar EUV 
(-50 to 1000 A) flux, which is a fundamental parameter to 
model physics, chemistry and dynamics of the upper atmo- 
sphere of planets. Since observations of solar EUV irradi- 
ance are not frequent and generally not available simultane- 
ously with the observation for the upper atmospheric stud- 
ies, the solar EUV flux model become important for mod- 
elling of aeronomical quantities in planetary atmospheres. 
Generally, solar EUV flux models are bin-averaged into 
numbers of wavelength bands and important solar emis- 
sion lines appropriate for the calculation of photoionization 
and photoelectron impact production rates. Characterisa- 
tion of the solar EUV flux for use in aeronomical and iono- 
spheric studies were developed during the seventies based 
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on the Atmospheric Explorcr-E (AE-E) data (Hintereggcr. 
1976; Hinteregger et al., 1981; Torr and Torr, 1985). Two 
AE-E reference spectra SC#21REF and F79050N have 
been pubhshed by Torr and Torr (1985) at 37 wavelength 
bins for solar minimum and maximum conditions, re- 
spectively. Later, based on the measured photoelec- 
tron flux, the short wavelength fluxes were increased by 
Richards et al. (1994), and they incorporated modified 
F74113 solar EUV flux in their EUVAC model. Detailed 
discussion on the development of solar EUV flux models is 
beyond the scope of this study; Lean (1990), Richards et al. 
(1994), and Lean et al. (2011) have provided reviews on so- 
lar EUV flux models. 

For a given solar activity there arc significant differences 
between the EUV fluxes reproduced by different solar flux 
models. These models are based on the different input pa- 
rameters and proxies, e.g., solar 10.7 cm radio flux (hence- 
forth referred to as solar F10.7) is used as the measure of so- 
lar activity, and used for paramctrization of solar EUV flux 
models (Richards et al., 1994; Tobiska and Barth, 1990). 
EUVAC model of Richards et al. (1994) is based on solar 
F10.7 and its 81-day average and also on the F74113 solar 
flux (in the EUVAC model, the F74113 flux below 250 A, 
and below 150 A, are doubled and tripled, respectively). 
SOLAR2000 model of Tobiska et al. (2000) is based on 
measured solar flux irradiance and various proxies and pro- 
vides solar flux from X-rays to infrared wavelengths, i.e., 
1-1000000 A. 

Different solar EUV flux models have been used to 
study the solar radiation interaction with upper atmo- 
sphere of Mars. Presently, SOLAR2000 (S2K) model of 
Tobiska (2004) and EUVAC model of Richards et al. (1994) 
are commonly used solar EUV flux models in aeronomical 
studies of Mars; e.g., Simon et al. (2009) and Huestis et al. 
(2010) have used EUVAC model, while Shematovich et al. 
(2008) and Cox et al. (2010) have used S2K model for the 
dayglow calculations on Mars. 

Fox et al. (1996) have studied the effect of different so- 
lar EUV flux models on calculated electron densities for low 
and high solar activity conditions. They have used 85315 
and 79050 solar fluxes of Tobiska (1991) and SC#21REF 
and F79050N AE-E reference solar spectra of Hinteregger 
(Torr et al., 1979) for low and high solar activity condi- 
tions, respectively. Fox et al. (1996) found that due to 
smaller fluxes at short wavelength range (18-200 A) m 



Hinteregger spectra, lower peak in electron density proflle 
is signiflcantly reduced (30-35%) compared to that calcu- 
lated using solar fluxes of Tobiska (1991). Buonsanto et al. 
(1995) calculated ionospheric electron density using EU- 
VAC model of Richards et al. (1994) and EUV94X solar 
flux model of Tobiska (1994). They found that photoioniza- 
tion rate in F2 region calculated by using EUV94X model 
is larger than that calculated using EUVAC model due to 
the large EUV fluxes in 300-1050 A wavelength range in 
EUV94X solar flux model. Recently, Jain and Bhardwaj 
(2011) have studied the effect of solar EUV flux models on 
N2 VK band intensities in Martian dayglow and showed 
that EUV flux models does affect the N2 dayglow emis- 
sions. Similar conclusion have been drawn for N2 dayglow 
emission on Venus (Bhardwaj and Jain, 2011) 

The aim of the present study is to calculate the impact 
of solar EUV flux models on CO Cameron band and COj 
UV doublet band intensities in Martian dayglow. We have 
used 37 bin EUVAC model of Richards et al. (1994) and 
S2K version 2.36 of Tobiska (2004) as the solar EUV fluxes. 
In these models, bins consist of band of 50 A width each and 
few prominent solar EUV lines, and are sufficient for the 
modelling of photoionization and photoelectron flux calcu- 
lations (Richards et al., 1994; Simon et al., 2009). Photo- 
electron flux, volume excitation rates and overhead inten- 
sities are calculated using both the solar EUV flux mod- 
els for low, moderate, and high solar activity conditions. 
Line of sight intensities for Cameron band and UV dou- 
blet emissions are calculated and compared with the latest 
observations by SPICAM instrument. 

2 Model 

Model atmosphere consists of five gases (CO2; N2, CO, O, 
and O2). Model atmosphere for solar minimum condition 
is taken from the Mars Thcrmospheric General Circula- 
tion Model (MTGCM) (Rougher et al., 1999, 2000, 2004). 
Model atmosphere for high solar activity is taken from Fox 
(2004). Photoabsorption and photoionization cross sec- 
tions for gases considered in the present study are taken 
from Schunk and Nagy (2000), and branching ratios for 
ionization in different states are taken from Avakyan ct al. 
(1998). 

Production mechanisms for CO(a"^H) are photon and 
electron impact dissociative excitation of CO2, electron dis- 
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sociativc recombination of COj , and electron impact ex- 
citation of CO. Since — > a'^H is a forbidden transi- 
tion, resonance fluorescence of CO is not an effective excita- 
tion mechanism. CO(a'^n) is a metastable state; Lawrence 
(1972) has measured its lifetime as 7.5 ± 1 ms, consis- 
tent with the measurement of Johnson (1972). Individual 
band lifetime can vary, e.g., lifetime of CO(a'^n) 1^ = 
level is around 3 ms (Gilijamse et al., 2007; Jongma et al., 
1997). Due to its long lifetime, cross section for the produc- 
tion of Cameron bands due to electron impact dissociation 
of CO2 is difficult to measure in the laboratory. Over- 
all, Cameron band cross section can have an uncertainty 
of a factor of 2 (Bhardwaj and Jain, 2009). Ajello (1971) 
reported relative magnitudes of the cross section for the 
(0, 1) band at 215.8 nm. Erdman and Zipf (1983) esti- 
mated the total Cameron band emission cross section of 
2.4 X 10""'^^ cm^ at 80 eV, whereas Ajello estimated a value 
of 7.1x10"" cm^ at 80 eV. Bhardwaj and Jain (2009) have 
fitted the e-C02 cross sections producing CO Cameron 
band based on the estimated value of Erdman and Zipf 
(1983). In our study cross section for Cameron band 
production due to electron impact on CO2 is taken from 
Bhardwaj and Jain (2009). Cross section for photodissoci- 
ation of CO2 producing CO(a^n) is taken from Lawrence 
(1972). To calculate the production rate of CO(a'^n) 
due to dissociative electron recombination process we 
have calculated the density of electron and major ions 
by including ion-neutral chemistry in the model. The 
coupled chemistry model is similar to that adopted in 
our other studies (Bhardwaj et al., 1996; Bhardwaj, 1999; 
Haider and Bhardwaj, 2005). Rate coefficients for ion- 
neutral reactions are taken from Fox and Sung (2001). Re- 
cently, Seiersen et al. (2003) have measured recombination 
rates for the e-COj collision, with yield of 0.87 for the 
channel producing CO of which CO(a'^n) production yield 
is 0.29 (Skrzypkowski et al., 1998; Rosati et al., 2003). Ion 
and electron temperatures are taken from Fox (2009). Ion 
and electron densities are calculated under steady state 
photochemical equilibrium. To calculate the density of 
CO(a'^n) we have also included various sources of loss of 
CO(a'^n) in our coupled chemistry model, which are given 
in Table 1. 

Major production sources of CO;^ in B state are pho- 
ton and electron impact ionization of CO2. Cross section 
for the formation of COj(B'^S]+) state due to electron im- 



pact ionization of CO2 is taken from Itikawa (2002), and 
cross section due to photoionization of CO2 is based on 
the branching ratio taken from Avakyan et al. (1998). For 
other gases electron impact cross sections have been taken 
from Jackman et al. (1977). 

Solar EUV flux has been taken at 1 AU (based on solar 
F10.7 at 1 AU as seen from the Mars, taking the Sun- 
Earth-Mars angle into consideration) and then scaled to 
the Sun-Mars distance. Fig. 1 shows the output of EUVAC 
and S2K solar EUV flux models for both solar minimum 
(20 July 1976) and solar maximum (2 August 1969) condi- 
tions at 1 AU. There are substantial differences in the solar 
EUV fluxes of EUVAC and S2K models; moreover, these 
differences are not similar in solar minimum and maximum 
conditions. In both, solar minimum and maximum condi- 
tions, solar flux in bands is higher in S2K than in EUVAC 
over the entire range of wavelengths, except for bins below 
250 A (150 A for solar minimum condition), whereas solar 
flux at lines is higher in EUVAC model for entire wave- 
length range. Overall higher solar fluxes above 250 A in 
S2K cause more photoionization. Higher photon fluxes be- 
low 250 A (d uring solar maximum condition) in EUVAC 
produce more higher energy electrons in the Martian atmo- 
sphere causing secondary ionizations (cf. Fig. 2) that can 
compensate for the higher photoionization in S2K model. 
One big difference between solar EUV fluxes of S2K and 
EUVAC models is the solar flux at bin containing 1026 
A H Ly-/3 line, which in both solar minimum and maxi- 
mum conditions is around an order of magnitude higher 
in S2K compare to EUVAC solar flux model (cf. Fig. 1). 
Flux at these wavelengths does not contribute to the pho- 
toionization, but are very important for dissociative exci- 
tation processes. Cross section for the photodissociation of 
CO2 producing Cameron band lies in longer (700-1080 A) 
wavelength regime (Lawrence, 1972). Hence, yield of pho- 
todissociation excitation of CO2 producing CO(a'^H) state 
would be larger when S2K solar EUV flux model is used. 

To calculate the photon degradation and generation of 
photoelectrons in the atmosphere of Mars, we have used 
Lambert-Beer law. Solar zenith angle (SZA) is 45° un- 
less mentioned otherwise in the text. Photoelectron en- 
ergy degradation and production rates of excitation of CO 
Cameron band and COj UV doublet band in the Martian 
atmosphere are calculated using Analytical Yield Spectrum 
(AYS) technique, which is based on the Monte Carlo model 
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(cf. Singhal and Bhardwaj, 1991; Bhardwaj and Singhal. 
1993; Bhardwaj and Michael, 1999a,b; Bhardwaj and Jain, 
2009). Details of calculation of photoelectron produc- 
tion rates and photoelectron flux have given in our 
earlier papers (Bhardwaj et al., 1990; Bhardwaj, 2003; 
Michael and Bhardwaj, 1997; Jain and Bhardwaj, 2011; 
Bhardwaj and Jain, 2011). 

Below 70 eV, photoelectron flux calculated using S2K 
is higher compared to that calculated using EUVAC model 
for low solar activity condition (Fig. 2). Above 70 eV, pho- 
toelectron flux calculated using EUVAC model is higher 
than that calculated using S2K model due to the larger 
solar EUV fluxes at shorter wavelength (< 200 A) in EU- 
VAC model (cf. Fig. 1). During solar maximum condition 
photoelectron fluxes calculated by using EUVAC model is 
higher than that calculated using S2K model (cf. Fig. 2) 
due to higher solar EUV flux at wavelength below 250 A 
in EUVAC model. During solar minimum, except at peaks 
and energies above 70 cV, photoelectron flux calculated 
using S2K model is around 1.4 times higher than that cal- 
culated using EUVAC model. During solar maximum con- 
dition photoelectron flux calculated using EUVAC model 
are higher than that calculated using S2K model. Photo- 
electron fluxes calculated using both solar EUV flux models 
are similar in shape but peaks around 20-30 eV are more 
prominent in electron flux calculated using EUVAC due to 
the higher solar EUV flux at lines (e.g.. He II Lyman-a 
line at 303.78 A) in EUVAC model (Fig. 1). The peaks in 
the 20-30 eV region of the photoelectron flux arise due to 
the ionization of CO2 in different ionization states by solar 
EUV flux at 303.78 line (Mantas and Hanson, 1979). Our 
calculated photoelectron fluxes are in good agreement with 
other model calculations (cf. Jain and Bhardwaj, 2011). 

Volume excitation rate is calculated for important pro- 
cesses producing CO(a'^n) and COj(B^I]+) states using 
photoelectron flux as 

V,{Z)=n{Z) cf,{Z,E)a,{E)dE, (1) 

where n{Z) is the density at altitude Z, <Ji{E) is the elec- 
tron impact cross section for i*'' process, for which thresh- 
old is Eth, and (f){Z,E) is the photoelectron flux. 



3 Results and discussion 

3.1 Low solar activity condition 

Wc run the model for low solar activity condition (similar 
to Viking landing), and calculated results arc compared 
with those of Fox and Dalgarno (1979) by taking the simi- 
lar model atmosphere. Model atmosphere is based on den- 
sity data derived from Viking 1 (Nier and McElroy, 1976). 
The Sun-Mars distance Dg_M = 1.64 AU and solar zenith 
angle is taken as 45°. 

Fig. 3 shows the calculated densities of COj and 
in the Martian upper atmosphere. The density of COj 
around peak and above calculated using S2K model is 
'-^30% higher than that calculated using EUVAC, which 
is due to higher production rate of COj when S2K model 
is used. Below 120 km, ion densities calculated by using 
EUVAC model are higher due to the higher photoelectron 
fluxes above 70 eV (cf. Fig. 2). There is a small disconti- 
nuity in the density of ion around 180 km, which is due 
to the sudden change in the electron temperature at 180 
km (see Fig. 2 of Fox, 2009). Our calculated ion densities 
are consistent with calculations of Fox (2004). 

Fig. 4 (upper panel) shows the proflles of production 
mechanisms of CO(a'^n) calculated using EUVAC and S2K 
solar EUV flux models. Around the peak of CO(a^n) pro- 
duction, the major source is photoelectron impact disso- 
ciation of CO2, while at higher altitudes photodissocia- 
tion excitation of CO2 takes over. Dissociative recombi- 
nation is about 13%, while photoelectron excitation of CO 
contribute about 3% to the total Cameron band excita- 
tion at the peak. The shape of volume excitation rate 
(VER) proflles and the altitude at the peak remain the 
same for all processes for the two solar flux models. How- 
ever, the magnitude of VERs calculated using S2K model 
arc about 40% higher than those calculated using EUVAC 
model. Contribution of electron impact dissociation of CO2 
producing CO(a^H) is higher in our studies than that in 
Fox and Dalgarno (1979). This is due to the higher value of 
CO(a'^H) production cross section in e-C02 collision used 
in the present study (having the value of 1 x 10^^^ cm^ at 
27 cV); Fox and Dalgarno (1979) used the cross section de- 
rived from Freund (1971) (having the value of 4 x 10^^^ cm^ 
at 27 cV). Due to larger photon flux at longer (700-1050 
A) wavelengths (region where photodissociation of CO2 
becomes important) in S2K model compared to EUVAC 
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model (of. Fig. 1), CO(a'^n) production due to photodis- 
sociative excitation is higher by ^50% when S2K model 
is used. Production rates of CO Cameron band for differ- 
ent processes calculated at the peak along with the peak 
altitude are given in Table 2. 

It is also clear from upper panel of Fig. 4 that the al- 
titude where the photodissociation of CO2 takes over elec- 
tron impact dissociation of CO2 in CO(a'^n) formation is 
slightly higher when EUVAC model is used (167 km for EU- 
VAC and 160 km for S2K solar flux model). In our model 
calculations, as well as in the work of Simon et al. (2009), 
photodissociation process becomes the major source at 
higher altitudes (> 160 km) and is a factor of 2 higher 
than the electron impact dissociation of CO2. 

For the COj UV doublet band, we have considered 
only photoionization and electron impact ionization of CO2 
producing COj in the B'^Sj state. Contribution of so- 
lar fluorescent scattering is very small, less than 10% (cf. 
Fox and Dalgarno, 1979), and hence it is not taken into ac- 
count in the present study. While calculating the emission 
from COj UV doublet, we have assumed 50% crossover 
from B to A state (Fox and Dalgarno, 1979). Fig. 4 (bot- 
tom panel) shows the production rates for COj in B state. 
Production of COj(B^E+) due to photoionization of CO2 
is about a factor of 3-4 higher than due to photoclcctron 
impact ionization. Here also we find that production rates 
calculated using S2K are higher than those calculated using 
EUVAC flux by about 50%, but peak altitude remains the 
same in both cases. Production rates calculated at peak 
along with peak altitude for COj UV doublet band are 
given in Table 2. 

Volume emission rates are height-integrated to calculate 
overhead intensities, which are presented in Table 3 for CO 
Cameron and UV doublet bands. For Cameron band, con- 
tribution of e-C02 process is maximum (64%) followed by 
photodissociation of CO2 (21%). Contributions of dissocia- 
tive recombination and e-CO process are around 10% and 
3%, respectively. For UV doublet band, the major contri- 
bution is coming from photoionization of CO2 (80%), the 
rest is due to electron impact ionization of CO2. 

To compare the model output with observed emissions 
we have integrated volume emission rates along the line of 



sight. Limb intensity at each tangent point is calculated as 

00 

/ = 2 y V(r)dr, (2) 


where r is abscissa along the horizontal line of sight, and 
V(r) is the volume emission rate (in cm^'^ s^^) at a par- 
ticular emission point r. The factor of 2 multiplication 
comes due to symmetry along the line of sight with respect 
to the tangent point. While calculating limb intensity we 
assumed that the emission rate is constant along local lon- 
gitude/latitude. For emissions considered in the present 
study the effect of absorption in the atmosphere is found 
to be negligible. 

Fig. 5 shows the calculated limb profiles of the CO 
Cameron and COj" UV doublet bands along with the 
SPICAM-observed limb intensity taken from Simon et al. 
(2009). Observed values are averaged over for the orbits 
close to Viking 1 condition (Ls'-^100-140°), low solar ac- 
tivity, and for SZA=45°. Below 100 km there is a sud- 
den increase in the observed intensity of CO J UV doublet 
band, which, according to Simon et al. (2009), is due to the 
very significant solar contamination below 100 km. Limb 
intensities calculated using S2K flux are '^40-50% higher 
than those calculated using EUVAC: clearly showing the 
effect of input solar EUV flux on the calculated intensi- 
ties. Magnitude of the calculated limb intensity of UV 
doublet band is in agreement with the observation, but for 
CO Cameron band calculated intensity is higher than the 
observed profile. For both bands, the calculated intensity 
profile peaks at higher ('^5 km) altitude in comparison with 
the observation-indicating a denser neutral atmosphere in 
our model. 

The dashed curves in Fig. 5 show intensities calculated 
after reducing the CO2 density by a factor of 1.5; a good 
agreement in the altitude of peak emission is seen between 
calculated and observed limb profiles. Though the reduc- 
tion in CO2 density shifted the altitude of peak emission 
downwards, the magnitude of calculated Cameron band in- 
tensity is still larger than the intensity measured by SPI- 
CAM. As pointed out in Section 2, the e-C02 cross sec- 
tions producing Cameron band are uncertain by a fac- 
tor of '^2. The calculated limb intensity profile for re- 
duced e-C02 cross section by a factor of 2 is also shown in 
Fig. 5. Cameron band intensities obtained after reducing 
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the density and cross section are in relatively close agree- 
ment with the observed values. In the model calculations 
of Simon et al. (2009) also the Cameron band intensity and 
its peak emission altitude were higher than the SPICAM 
observed values. They have to reduced the density of CO2 
by a factor of 3 and e-C02 Cameron production cross sec- 
tion by a factor of 2 to bring their calculated intensity 
profile in agreement with the SPICAM observation. 

3.2 SPICAM observations 

Leblanc et al. (2006) have presented detailed analysis of 
SPICAM data during the period October 2004 to March 
2005, spanning the solar longitude (Ls) from 101° to 171°. 
They divided the total data set in two periods of so- 
lar longitude: first, Ls = 101° to 130°, and second, Ls 
= 139° to 171° (cf. Table 2 of Leblanc et al., 2006). 
Leblanc et al. found that the altitude of peak emission for 
CO2 UV doublet and CO Cameron bands is around 10 
km higher for Ls > 138° (122.5 km and 132.5 km for UV 
doublet and Cameron bands, respectively) compared to Ls 
< 130° (112.5 km and 117.5 km, for the same emissions). 
Leblanc et al. (2006) could not provide the reason for the 
higher altitude of peak emission for Ls > 130° observa- 
tions. Later, Forget et al. (2009) derived neutral densities 
in Martian upper atmosphere using the SPICAM instru- 
ment in stellar occultation mode for the same observation 
period. Forget et al. (2009) found that their is a sudden 
increase in the CO2 density in the Martian upper atmo- 
sphere for Ls ^ 130°-140°, which they attributed to a dust 
storm. Dust storm can heat the lower atmosphere and thus 
increase the densities at higher altitudes, which could ex- 
plain the higher altitude for peak emission observed by the 
SPICAM for Ls > 130° observations. Increase in the alti- 
tude of peak intensity of dayglow emissions clearly shows 
the effect of dust storms on Martian dayglow emissions. 

Comparisons of SPICAM observations with model 
calculated dayglow emissions have been performed 
by Shematovich et al. (2008), Simon et al. (2009), and 
Cox et al. (2010). Simon et al. (2009) have used one 
dimensional Trans-Mars model with EUVAC solar flux 
model, whereas Shematovich et al. (2008) and Cox et al. 
(2010) have used Monte Carlo model with S2K solar fiux. 
In the present study, we have taken both EUVAC and S2K 
models and have calculated the Cameron band and UV 
doublet band emissions using the Analytical Yield Spec- 



trum method; the results are compared with the SPICAM 
observations. 

3.2.1 First Case (Ls < 130°) 

To model the SPICAM observations for Ls < 130° the 
model atmosphere is based on MTGCM of Bougher et al. 
(1999) (taken from Shematovich et al., 2008). Calculations 
are made for MEX orbit no. 983 (24 Oct. 2004) when 
Ds-M = 1.64 AU, F10.7 = 87.7 (F10.7A = 107.3). 

Fig. 6 (upper panel) shows the volume excitation rate 
of C0(a3n). The total VER calculated using S2K fiux is 
only slightly higher than that obtained using EUVAC flux. 
However, the production rate due to photodissociative ex- 
citation of CO2 is around 50% higher when S2K model 
is used. Another interesting feature is the dissociative re- 
combination (DR) process, whose contribution is '^18% in 
the total intensity, roughly equal to the photodissociative 
excitation process (DR contribution is even higher than 
photodissociative excitation around production peak when 
EUVAC model is used), and it is signiflcantly higher than 
compared to the DR process in Viking condition case (see 
Table 3). This is due to the higher density of CO J ion 
compared to Viking condition (see Fig. 3). Leblanc et al. 
(2006) mentioned that higher values of CO J can contribute 
up to 30% to the Cameron band production depending on 
the solar zenith angle. To account for DR in Cameron 
band production, Shematovich et al. (2008) and Cox ct al. 
(2010) have taken COj and electron densities from Fox 
(2004) for low solar activity condition. Since SPICAM ob- 
servations are made during moderate solar activity condi- 
tion, the contribution of DR in Cameron band production 
would be lower in their calculations. 

Fig. 6 (bottom panel) shows the production rates of 
CO2 UV doublet band. Total rate calculated using both 
solar flux models is peaking at same altitude (^125 km), 
but total production rate calculated using S2K model is 
higher (around 10%) than that calculated using EUVAC 
model. 

Table 3 shows the overhead intensities of CO Cameron 
and CO2 UV doublet bands calculated using both EU- 
VAC and S2K solar fiux models. Contribution of different 
processes in Cameron band production is slightly different 
than that in the low solar activity condition. Contribution 
from CO2 photodissociation is slightly reduced (17%, 21% 
when S2K model is used) , while dissociative recombination 
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contribution is increased ('^18%). Contribution of e-C02 
and e-CO processes remains almost same; 61 (64%) and 4 
(4%), respectively, when calculated using EUVAC (S2K) 
model. For UV doublet band photoionization of CO2 re- 
mains the dominant process contributing around 80% to 
the total overhead intensity. 

We have also calculated overhead intensities of major 
vibrational bands of Cameron system, which have been 
observed in Martian dayglow, using Frank-Condon fac- 
tors from Halmann et al. (1966) and branching ratio from 
Conway (1981). Table 4 shows the calculated overhead in- 
tensities of major vibrational bands of CO Cameron band 
system. Contribution of major vibrational bands to the 
total overhead Cameron band intensity is around 10, 10, 
16, and 8% for (0, 0), (0, 1), (1, 0), and (2, 0) bands, 
respectively. 

Fig. 7 shows the limb intensity profiles of Cameron 
and CO2 UV doublet bands. SPICAM-observed intensi- 
ties of Cameron and UV doublet bands averaged over Ls — 
100-130° observations (Leblanc et al., 2006) are also shown 
in the Fig. 7. Limb intensities of COj UV doublet and 
Cameron bands calculated by using S2K model are ^6% 
and ~15%, respectively, higher compared to those obtained 
using EUVAC model. Calculated intensities for both solar 
flux models are higher than the SPICAM-observed values. 
In analogy to the Viking case, we reduced the e-C02 cross 
sections producing Cameron band. The resulting intensity 
profile (also shown in Fig. 7) is still higher than the obser- 
vation around emission peak. Calculated intensity of CO J 
UV doublet band is also higher near the peak emission than 
the observed intensity. Altitude of the calculated intensity 
for both CO Cameron and UV doublet bands peaks ^2 
to 3 km higher than the observations, which is well within 
the observational uncertainties. Line of sight intensity of 
different vibrational transitions of Cameron band at the 
altitude of peak emission, which is ~120 km for first case, 
are shown in Table 4. 

Fig. 8 shows the calculated intensity ratio of UV dou- 
blet to Cameron band along with the observed ratio de- 
rived from SPICAM observations (Leblanc et al., 2006). 
At lower altitudes calculated ratio is in agreement with 
observation (^0.18). The ratio remains almost constant 
up to ^120 km (where Cameron band and UV doublet 
emission peaks), starts gradually decreasing with altitude 
and becomes almost constant after 150 km. The observed 



ratio decreases almost monotonically from 100 km all the 
way to 180 km. Leblanc et al. (2006) have not found any 
dependence of SZA on the UV doublet to Cameron band 
intensity ratios, though they have observed a weak depen- 
dence of this intensity ratio on the solar activity. From 
the observed intensity ratio profile it is clear that in upper 
atmosphere Cameron band intensity is increasing steadily 
compare to UV doublet band, which indicates a difference 
in the source of production of Cameron band and UV dou- 
blet band. That source could be the dissociative recom- 
bination process which is sensitive to the density of CO^ 
ion (as shown in the Fig. 3). Loss of CO J ions at higher 
altitudes (> 200) can reduce the intensity of UV doublet 
and hence decreases the intensity ratio value. 

3.2.2 Second Case (Ls > 130°) 

As discussed earlier, due to dust storm during SPICAM ob- 
servations for Ls greater than 130°, atmospheric densities 
were higher resulting in altitude of peak emission shifting 
to higher altitudes (^132.5 km for Cameron band emis- 
sion). For Mariner 6 and 7 observations the intensity of 
CO Cameron band peaked at altitude of ~133 km. Mariner 
observations were carried out during solar maximum con- 
dition (F10.7 ^ 180), whereas SPICAM observations are 
made during moderate solar activity condition. To model 
dayglow emissions for Ls > 130°, we have made calcula- 
tion for MEX orbit 1426 (26 Feb. 2005), taking model 
atmosphere from Fox (2004) for high solar activity condi- 
tion. Sun-Mars distance is 1.5 AU, F10.7 = 98 (F10.7A 
= 97). The EUV flux at 1 AU calculated using EUVAC 
model remains the same for first (Ls < 130°) and second 
(Ls >130°) cases. This is because in the EUVAC model 
average of F10.7 and F10.7A (81-day average) is used to 
scale the solar flux, and on both days average of F10.7 
and F10.7A does not change (it is 97.5 on both days), al- 
though the F10.7 flux increased by 10 unit. S2K model 
does not depend on the F10.7 alone, but on other prox- 
ies also (Tobiska, 2004), hence flux calculated using S2K 
model is different on the two days. 

Fig. 9 (upper panel) shows the VER of CO(a^n) cal- 
culated using S2K and EUVAC models. Total VER calcu- 
lated using S2K is about 17% higher than that calculated 
using EUVAC model. Major differences are in the CO(a'^n) 
production due to dissociative excitation of CO2 by pho- 
ton and electron impact, which are more than 50% and 10% 
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higher, respectively, when S2K model is used. Total pro- 
duction rate of Cameron band (calculated using EUVAC 
model) maximises at an altitude of 134 km with a value of 
about 3328 cm"'^ s"-'^, which is around 10 km higher than 
that in the first case (Ls < 130°). Although production 
rate (3528 cm~'^ s~^) at peak altitude is higher in the first 
case, but at higher altitudes rate increases faster in the 
second (Ls > 130°) case, e.g., at 200 km, Cameron band 
production rate is 61 em^'^ s^^ in second case, whereas in 
first case it is only 3 cm^^ s^^. In both, first and second 
cases, for EUVAC model, the altitude where photodisso- 
ciation of CO2 takes over electron impact dissociation of 
CO2 is around 30% higher than that for S2K model. 

Bottom panel of Fig. 9 shows production rates of CO;^ 
UV doublet band. Total excitation rate calculated using 
the S2K model is about 12% higher than that calculated 
using EUVAC model. Similar to the Cameron band, COj 
UV doublet production rate at peak is lower than that 
in the first case, but at higher altitudes UV doublet pro- 
duction rate becomes higher in the second case. Table 3 
shows the overhead intensities of Cameron and COj UV 
doublet bands. Contribution of photodissociation of CO2; 
e-C02, DR, and e-CO processes to the total Cameron band 
production is 16 (22%), 62 (57%), 11 (11%), and 9 (8%), 
respectively, when EUVAC (S2K) model is used. 

Fig. 10 shows the calculated limb intensity of UV dou- 
blet and Cameron band along with SPICAM-observed in- 
tensities of Cameron band and CO J doublet band for MEX 
orbit 1426 on 26 Feb. 2005 (Shematovich et al., 2008). In- 
tensities calculated using S2K model are higher by ^12- 
18% than those calculated using EUVAC model. Altitude 
of peak emission of calculated and observed intensity pro- 
files is in good agreement with each other (e.g., ^128 km 
for Cameron band) within the uncertainties of observations 
and model calculations. However, intensities calculated us- 
ing both solar flux models are higher than the observations. 
Calculated intensities of CO Cameron band after reducing 
the e-C02 cross section by a factor of 2 are also shown 
in Fig. 10, which is in good agreement with the observed 
values. However, the calculated intensity of UV doublet 
band is 20-30% higher than the SPICAM-observed val- 
ues. Uncertainties in photoionization cross sections for the 
production of COj(B^E+) can also be one of reasons for 
the differences between observed and calculated intensity 
of C0+(B2I]+) UV doublet emission. 



In our calculation (using EUVAC model and reduced 
e-C02 cross section for Cameron band production), peak 
brightness of Cameron band is 249 (255) kR with peak lo- 
cated at 120 (128) km in first (second) case. For COj UV 
doublet band, peak brightness is 31 (33) kR with peak lo- 
cated at 119 (127) km in first (second) case. Altitude of 
peak brightness remains the same when S2K model is used 
but intensities of Cameron and COj UV doublet bands 
for first (second) case are about a factor of 1.1 (1.18) and 
1.07 (1.11), respectively, higher than that calculated using 
EUVAC model. 

3.3 Solar Maximum (Mariner Observa- 
tions) 

First dayglow measurements on Mars were carried out by 
Mariner series of spacecraft. Mariner 6 and 7 observations 
were taken during the solar maximum conditions (July- 
August, 1969; F10.7 = 186 at 1 AU), and are the only 
Martian dayglow observations so far taken during the so- 
lar maximum condition. Mars was also at perihelion (dis- 
tance between Sun and Mars was around 1.42 AU) during 
Mariner observations. We run our model for the condition 
similar to the Mariner observation. Model atmosphere for 
solar maximum condition is taken from Fox (2004). 

Fig. 1 (bottom panel) shows the solar EUV flux during 
Mariner observations (2 Aug. 1969, F10.7 = 180) from EU- 
VAC and S2K models. As discussed in Section 2, overall 
the solar EUV flux is higher in S2K model, except at wave- 
lengths below 250 A where EUVAC flux is larger. Fig. 2 
(bottom panel) shows photoelectron fluxes calculated at 
different heights in solar maximum condition using both 
solar flux models. Unlike in the solar minimum condition, 
the ratio of photoelectron flux calculated using S2K to that 
of EUVAC is less than 1 at most of the energies. At elec- 
tron energies below 30 eV, fluxes calculated using S2K and 
EUVAC models are almost equal at altitudes of 130 and 
160 km. Above 30 eV, photoelectron fluxes calculated us- 
ing EUVAC model are higher than those calculated using 
S2K model. The higher photon fluxes in EUVAC model 
at shorter wavelengths below 250 A [cf. Fig. 1] produce 
higher energy photoelectrons that can produce more sec- 
ondary electron through ionization, and hence compensates 
for the higher solar EUV flux above 250 A in the S2K 
model. Table 3 shows the calculated overhead intensities 
of Cameron band and UV doublet band using EUVAC and 
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S2K models. Cameron band production due to electron 
impact dissociation of CO2 is higher when EUVAC model 
is used, which is due to the higher photoelectron fluxes. 
Contribution of photodissociation excitation calculated us- 
ing S2K model is higher, due to higher EUV flux at longer 
wavelengths (specially flux in the 1000-1050 A bin). 

Fig. 11 (upper panel) shows the VER of CO(a'^n) for 
higher solar activity condition, calculated using EUVAC 
and S2K solar flux models. Due to the higher photoelec- 
tron flux, CO(a^n) production due to e-C02, e-CO, and 
dissociative recombination are higher when EUVAC model 
is used. Photodissociative excitation of CO2 producing 
Cameron band for S2K model is still higher by ~50% than 
for EUVAC model, which is due to the higher EUV fluxes 
at longer wavelengths in the S2K model. Similar to that in 
the previous cases, the cross over point between photodis- 
sociation and electron impact dissociation of CO2 forming 
CO(a'^n) occurs at higher altitude when EUVAC model is 
used. Bottom panel of Fig. 11 shows the production rates 
of CO2 UV doublet band. Here also calculated values us- 
ing EUVAC model is slightly higher than that calculated 
using S2K model. 

During solar minimum condition total volume produc- 
tion rate of Cameron and UV doublet bands calculated us- 
ing S2K model is higher than that calculated using EUVAC 
model, whereas in solar maximum it is vice- versa. Except 
photodissociation excitation process producing Cameron 
band, production rates due to other processes calculated 
by using EUVAC model are higher than that calculated by 
using S2K model. In both, solar minimum and maximum 
conditions, Cameron band production due to photodisso- 
ciative excitation is about 50% higher, when S2K model is 
used. 

Fig. 12 shows model intensities of Cameron band and 
CO^ (B - X) ultraviolet doublet band calculated using both 
EUVAC and S2K models at SZA = 45° along with inten- 
sities observed by Mariner 6 and 7. Limb intensities mea- 
sured by Mariner 6 and 7 on Mars are at SZA = 27° and 
0°, and at SZA = 44° and 0°, respectively. Calculated 
limb intensities using EUVAC model at SZA = 0° are also 
shown in the Fig. 12. Limb intensities calculated using EU- 
VAC model are only slightly higher than those calculated 
using S2K model. There were no observations at emis- 
sion peak for both Cameron band and COj(B - X) band. 
Below the emission peak, there are few observations for 



Cameron band, but ultraviolet doublet observations were 
not available. Solar zenith angle effect is clearly visible at 
altitudes below 150 km, where intensity is larger and emis- 
sion peak shift deeper in the atmosphere for lower SZA. 
Calculated intensities of Cameron and UV doublet bands 
are lower than the observed values. Unlike previous cases, 
calculated intensities of Cameron and UV doublet band 
emissions at the altitude of peak emission are slightly lower 
than the observation. Stewart et al. (1972) has pointed out 
that due to calibration problem in Mariner 6 and 7 instru- 
ment the observed values can be higher. As in the previous 
cases a reduction in e-C02 cross section is required to get 
an agreement between observed and calculated intensity. 
Calculated intensity of Cameron band after reducing the 
e-C02 cross section by a factor of 2 is also shown in Fig- 
ure 12 for SZA=0°. 

3.4 CO(a3n) density 

Density of CO(a'^n) is calculated under photochemical 
equilibrium condition. Radiative decay is the dominant 
loss process of CO(a'^n), the contribution from other pro- 
cesses are negligible (cf. Bhardwaj and Raghuram, 2011). 
Fig. 3 shows the density of CO(a'^n) calculated using EU- 
VAC and S2K EUV flux models for low solar activity condi- 
tion. The calculated column density of CO(a"^n) molecule 
is 4.6 X 10^ (6.5 X 10^) cm^^ for the solar minimum condi- 
tion using EUVAC (S2K) solar EUV flux model. Except in 
the solar maximum condition, density of CO(a'^n) molecule 
calculated using S2K model is higher than that calculated 
using EUVAC model. During solar maximum condition, 
CO(a'^n) density calculated using EUVAC model is slightly 
higher at peak (around 5%), but at altitudes above 140 km, 
density calculated using S2K model becomes higher ('^10% 
at 200 km). 

The shape of the density of CO(a'^n) is similar to that 
of its production rate (cf. Fig. 4) since the main loss mech- 
anism of CO(a^n) is radiative decay whose value is inde- 
pendent of altitude. Hence, the density of CO(a'^n) in the 
Martian atmosphere can be represented by 

[C0(a3n)]^[^Q-](^^+^-) + [^Q^][-°]^- (3) 

K4 

where Ki, K2, K3, and K4 are as described in Table 1. 
Ki and K2 are photodissociation rate and electron impact 
dissociation rate of CO2, respectively, K3 is dissociative 
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recombination, is radiative dceay loss, and ne is the 
electron density. The values of Ki (photodissociation rate) 
in units of s^^ at the top of atmosphere in case of EU- 
VAC (S2K) model are 7.5 x 10"^ (1.1 x lO"'^), 8.7 x 10^^ 
(1.3 X 10-^), 1.03 X 10^^ (1.6 X 10-'^), and 1.5 x 10"^ 
(2.25 X 10~^) in the solar minimum, first case, second case, 
and solar maximum, respectively. 

4 Summary and Conclusion 

Present study deals with the model calculations of CO 
Cameron band and COj doublet ultraviolet emissions in 
Martian dayglow and the impact of solar EUV flux on 
the calculated intensities. Photoelectrons generated due to 
photoionization in the Martian atmosphere have been de- 
graded in the atmosphere using Monte Carlo model-based 
Analytical Yield Spectrum technique. Emission rates of 
Cameron and COj UV doublet bands due to photon and 
electron impact on CO2 have been calculated using EUVAC 
and S2K solar EUV flux models. Densities of prominent 
ions and CO(a^n) in Martian upper atmosphere are calcu- 
lated under steady state photochemical equilibrium condi- 
tion. Production rates of Cameron and COj UV doublet 
bands are height-integrated to calculate overhead intensity 
and along the line of sight to obtain limb intensities. Limb 
intensities are compared with the SPICAM/Mars Express 
and UV spectrometer/Mariner observed intensities. 

Due to higher EUV fluxes at longer (700-1050 A) wave- 
lengths in the S2K model, the contribution of photodisso- 
ciation of CO2 in producing Cameron band is about 50% 
higher in low as well as in high solar activity conditions. 
Variations in EUV fluxes at longer wavelengths from solar 
minimum to solar maximum are less prominent in the EU- 
VAC solar EUV flux model compared to the S2K model. 

For low solar activity condition, limb intensities of 
Cameron and CO^ UV doublet bands around peak bright- 
ness calculated using S2K model are ~3G-40% higher than 
those calculated using EUVAC model. Comparison of 
calculated intensities has been made with the SPICAM- 
observed values for condition similar to the Viking. In- 
tensities calculated using both S2K and EUVAC models 
are higher than the observed values. Calculated altitude of 
emission peak of CO Cameron and CO^ UV doublet bands 
is also higher by ~5 km than the observed value. A reduc- 
tion in the e-C02 cross section forming Cameron band by 



a factor of 2 and the density of CO2 in model atmosphere 
by a factor of 1.5 brings the calculated intensity (using EU- 
VAC model) of Cameron band in close agreement with the 
SPICAM observation. 

While modelling the recent observations made by SPI- 
CAM on-board Mars Express, we have taken two set of 
conditions with different model atmospheres and solar lon- 
gitudes. In the first case, Ls < 130° (Ls = 100-130°), at- 
mosphere is taken from Mars Thermospheric General Cir- 
culation model (Bougher et al., 1999, 2000, 2004) and cal- 
culations are made for the day 24 October 2004 with mod- 
erated solar activity flux (F10.7 = 88). Total intensities of 
CO Cameron and COj UV doublet bands calculated us- 
ing S2K model are around ^-^6-15% higher than those cal- 
culated using EUVAC model. Contribution of CO2 pho- 
todissociative excitation in Cameron band production is 
50% higher when S2K model is used. Dissociative recom- 
bination of C02^ is an important source of Cameron band 
in this case (cf. Fig. 6) due to higher densities of COj ion 
(Fig. 3) compared to those calculated for low solar activity 
condition (Viking condition). 

Calculated intensities of Cameron and UV doublet 
bands have been compared with the SPICAM-observed 
limb intensities. Intensities calculated using S2K and EU- 
VAC solar flux models are higher than the observed values 
by a factor of 1.7 to 2 for Cameron and a factor of 1.4 for 
UV doublet bands (see Fig. 7). We found that altitude of 
peak emission of both Cameron and UV doublet bands are 
2 to 3 km higher than observed profiles. This discrepancy 
in observed and calculated intensities and altitudes could 
be due to the fact that observed values are averaged over 
several days of observations while calculation are carried 
out for a particular day. 

Due to the dust storm during Ls > 130°, observed emis- 
sion peak is around 10 km higher for both Cameron and 
UV doublet compared with the SPICAM-observed values 
for MEX orbit 1426 on 26 Feb. 2005 (Shcmatovich et al., 
2008). To model the emission during Ls > 130°, we have 
taken atmospheric model for solar maximum condition. In- 
tensities calculated using the S2K solar flux model are ~8- 
18% higher than those calculated using the EUVAC model. 
These calculated intensities are higher than the observed- 
averaged values by a factor of ~2 for Cameron band and 
^50% for UV doublet band. The calculated intensity of 
Cameron band (after reducing the e-C02 cross sections by 
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a factor of 2) is in agreement with the observed values 
(Fig. 10). 

In all three conditions discussed above, i.e., low solar 
activity (Viking), and first (Ls < 130°) and second (Ls 
> 130°) cases, calculated intensities of both Cameron and 
UV doublet bands are higher than observations. On an av- 
erage, model values of Cox et al. (2010) for Cameron and 
CO J UV doublet bands are around a factor of 1.74 and 
1.41, respectively, higher than the SPICAM observations. 
Simon et al. (2009) also found that their calculated inten- 
sities of Cameron and UV doublet bands are around 25% 
higher than the SPICAM-observed values. This shows that 
these discrepancies in the model and observed values are 
due the uncertainties in the input physical parameters in 
the model. Uncertainties in cross sections, namely, e-C02 
cross section producing Cameron band and photoionization 
of CO2 forming UV doublet band can be one of the causes 
of discrepancies in the model and observations. 

Calculations are also made for the high solar (Mariner 
6 and 7 observations) activity condition. Though the con- 
tribution of CO2 photodissociative excitation in Cameron 
band production is higher when S2K model is used, but the 
total intensity of CO Cameron band calculated using the 
EUVAC model is slightly higher than that calculated using 
the S2K model. This is because of the higher photoelectron 
flux when EUVAC model is used (Fig. 2). The calculated 
intensities of both Cameron and UV doublet band are lower 
than the observed values (Fig. 12). 

Following conclusion can be drawn from the present 
study: 

1. Generally, solar EUV fluxes in bands are higher in 
S2K model except at few bands at shorter wavelength 
range (< 250 A). Solar EUV fluxes at Ion ger wave- 
lengths are higher in S2K model, specially in the 
1000-1050 A bin, where flux is around an order of 
magnitude higher than the corresponding flux in EU- 
VAC model. Solar EUV flux at lines is smaller in the 
S2K model compared to that in the EUVAC model. 

2. Due to higher EUV flux at hues in the EUVAC model, 
peaks at 20-30 eV range in the photoelectron flux are 
more prominent when EUVAC model is used. 

3. During high solar activity condition, calculated pho- 
toelectron fluxes are higher for EUVAC model due 
to higher EUV fluxes below 250 A in the EUVAC 



model. Hence, intensities calculated using EUVAC 
model are higher by 5-10% than those calculated us- 
ing S2K model. 

4. During solar minimum condition, the Cameron and 
CO2 UV doublet intensities calculated using S2K so- 
lar flux model are ~30-40% higher than those calcu- 
lated using the EUVAC model. 

5. During both, solar minimum as well as maximum 
conditions, Cameron band production due to pho- 
todissociative excitation of CO2 is about 50% higher 
when S2K solar EUV flux model is used. 

6. Altitude of peak production rate of Cameron and 
CO^ UV doublet bands is independent of solar EUV 
flux model used in the calculations. However, for 
the Cameron band the altitude where photodissoci- 
ation of CO2 takes over electron impact dissociation 
is higher in the EUVAC model compared to that in 
the S2K model. 

7. Reduction in the e-C02 cross section producing 
Cameron band and photoionization cross section pro- 
ducing CO^ UV doublet band is required to bring 
the model calculations in agreement with the obser- 
vations. 

8. For a given set of observation, and accounting for the 
uncertainties in the cross sections, intensities calcu- 
lated using the EUVAC model are in better agree- 
ment with the observation than those calculated us- 
ing the S2K model. 

Simultaneous observation of solar EUV flux with dayglow 
measurements would be very helpful in improving our un- 
derstanding of the processes that governs the dayglow emis- 
sions on Mars. More accurate measurements of cross sec- 
tions for electron impact dissociation of CO2 producing 
Cameron band and photoionization of CO2 in B state are 
required for the better modelling of CO Cameron band and 
CO^ UV doublet band in Martian atmosphere, as well as in 
other C02-containing atmospheres, like Venus and comets. 
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Table 1: Major reactions for the production and loss of CO(a^n). 



Reaction Rate (cm"^ s ^ or s ^) Reference 



CO2 + 


C0(a3n) + 


oc^p) 


Model (Ki) 


Lawrence (1972) 


CO2 + 


^ C0(a3n) -F 


- + c- 


Model (K2) 


Present work 


CO + e;, - 

C0+ + e- 


^ C0(a3n) + 


c~ 


Model 


Present work 


^ CO(a^^n) 4 







Seiersen et al. (2003); 
Rosati et al. (2003) 


C0(a3n) + 


CO2 ^ CO ^ 


- C02 


1.0 X 10-" 


Skrzypkowski et al. 
(1998) 


C0(a3n) + 


CO ^ CO + 


CO 


5.7 X 10-" 


Wysong (2000) 


C0(a3n) 


> CO + Iw 




K4= 1.26 X 10^ 


Lawrence (1971) 



%h ~ photoelectron. 

tKa = 6.5 X 10"^ (300/Te)°-* x 0.87 x 0.29 cm^ s"^; here 0.87 is yield of dissociative recombination of COj producing CO, and 
0.29 is yield of CO(a^n) produced from CO. 
Ki, K2, Ks, and K4, are described in equation 3. 



Table 2: Production rates (in cm ^ s ^) with peak altitude for low solar activity condition (similar to Viking condition). 
Production Process Cameron band COj UV doublet 

391 at 131 km 
95 at 127 km 



478 at 131 km 
%h — photoelectron. 



CO2 + hi/ 465 at 137 km 

CO2 + e^"^ 1470 at 128 km 

81 at 130 km 



ph 

Total 2094 at 131 km 



COJ + e- 267 at 133 km 
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Tabic 3: Overhead intensities (in kR) of CO Cameron and COj UV doublet bands. 



Emissions 



EUVAC 



SOLAR2000 



h!/+C02 C+CO2 DR9 e+CO Total 



h!/+C02 e+C02 DR e+CO Total 



Cameron band 
UV doublet 

Cameron band 
UV doublet 



Cameron band 
UV doublet 

Cameron band 
UV doublet 



1.2 
1.1 



1.4 
1.4 



1.6 
1.7 

2.3 
3.1 



3.7 
0.2 



5.4 
0.4 



6.1 
0.4 

11.8 
0.8 



Low solar activity (Viking condition) 

0.6 0.2 5.8 1.8 5 

1.3 1.4 0.3 

SPICAM/Mars-Express, First Case (Ls < 130°) 

1.6 0.4 8.8 2.1 5.7 

1.8 1.6 0.4 

SPICAM/Mars-Express, Second Case (Ls > 130°) 

1.1 0.9 9.8 2.6 6.7 

2.1 1.9 0.4 

Higher solar activity (Mariner observations) 

2.4 1.7 18.3 3.6 10.2 

3.9 2.8 0.7 



1.3 



2.2 



0.3 



0.4 



1.6 



8.2 
1.7 



10.1 

2 



11.7 
2.3 

17.7 
3.5 



s Dissociative recombination (e + CO2 



Table 4: Overhead and limb intensities of different vibrational bands of CO Cameron band system calculated after reducing 
the e-C02 cross section by a factor of 2 in the first case (Ls < 130°). 



Band 


Band Origin 


Overhead (R) 


Limb (kR)t 


{u' - ly") 


A 


EUVAC 


S2K 


EUVAC S2K 


- 


2063.0 


716 


845 


20 24 


- 1 


2158.4 


713 


840 


20 23.2 


- 2 


2261.0 


311 


367 


9 10 


- 3 


2374.0 


79 


93 


2 2.6 


1 - 


1992.5 


1151 


1358 


32 38 


2 - 


1927.5 


603 


711 


17 20 


3 - 


1868.1 


181 


213 


5 6 


4 - 


1813.0 


38 


45 


1 1.3 



Limb intensities are given at altitude (120 km) of peak emission. 



18 



Jain and Bhardwaj., 2011, PSS, doi:10.1016/j.pss.2011.08.010 



(a) 



S2K Solar Minimum 



28 

I 1 

X X 

X X 



iL 



EUVAC Solar Minimum 



50 150 250 350 



450 550 650 
Wavelength (A) 



750 850 950 1050 




50 150 250 350 450 550 650 750 850 950 1050 

Wavelength (A) 



Figure 1: Solar photon flux, in bands and at lines, as a funetion of wavelength in EUVAC and S2K solar EUV flux models, 
(a) for the low solar activity condition on July 1976 (similar to Viking landing, F10.7 — 70), and (b) for high solar activity 
condition on August 1969 (similar to Mariner 6 and 7 observations period, F10.7 = 186). 
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Figure 2: Calculated photoelectron fluxes for low (upper panel) and high (lower panel) solar activity conditions at SZA = 
45°. The ratio of the photoelectron flux at 130 km calculated using the two solar flux models is also shown with magnitude 
on right side Y-axis. Thin dotted horizontal line depicts the S2K/EUVAC ratio = 1. 
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Figure 3: Densities of COj and ions and CO(a"^n) molecule for solar minimum condition calculated using EUVAC 
(solid curve) and S2K (dashed curve) solar EUV flux models. Density of CO(a'^n) molecule is plotted after multiplying 
by a factor of 100. Dotted curves show the densities of COj and ions for first case (Ls < 130°) using EUVAC solar 
flux. 
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Figure 4: Calculated produetion rates of the CO(a'^n) (upper panel) and COj(B^I]+) (bottom panel) for low solar activity 
condition (Ls~100-140°). DR stands for dissociative recombination. 
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Figure 5: Calculated limb profiles of COj UV doublet bands (left panel) and CO Cameron (right panel) for low solar 
activity condition. Solid squares with error bars represents the SPICAM-observed values taken from Simon et al. (2009). 
Dashed curves show the calculated intensity (using EUVAC model) after reducing the density of CO2 by a factor of 1.5. 
Dash-dotted curve shows the calculated intensity (using EUVAC) of Cameron band with reduced density (by a factor of 
1.5) and reduced (by a factor of 2) e-C02 Cameron band production cross section. 
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Figure 6: Calculated production rates of the CO(a'^n) (upper panel) and COj(B^I]+) (bottom panel) for solar longitude 
Ls < 130°. DR stands for dissociative recombination. 
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Figure 7: Calculated limb profile of COj" UV doublet band (left panel) and CO Cameron band (right panel) for Ls < 
130°. Symbols represent the SPICAM-observed values taken from Leblanc et al. (2006). Dash-dotted curve shows the 
calculated intensity of Cameron band with reduced (by a factor of 2) e-C02 cross section. 
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Figure 8: Altitude variation of intensity ratio of UV doublet and CO Cameron band system. Calculated ratio is shown 
for EUVAC solar flux model with Cameron band production cross section in e-C02 collision is reduced by a factor of 2. 
SPICAM-observed ratio is from Fig. 9(a) of Leblanc et al. (2006). 
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Figure 9: Calculated production rates of the CO(a'^n) (upper panel) and COj(B^I]+) (bottom panel) for solar longitude 
Ls > 130°. DR stands for dissociative recombination. 
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Figure 10: Calculated limb profiles of CO^ UV doublet band (left panel) and CO Cameron band (right panel) for Ls 
> 130°. Dash dotted curve shows the Cameron band intensity for EUVAC model with e-C02 cross section forming 
CO(a^n) reduced by a factor of 2. Open circles with error bars represent the SPICAM-observed intensity taken from 
Shematovich et al. (2008). 
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Figure 11: Calculated production rates of the CO(a'^n) (upper panel) and COj(B^S]+) (bottom panel) for solar maximum 
condition. DR stands for dissociative recombination. 
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Figure 12: Calculated limb intensity for the higher solar activity condition similar to Mariner 6 and 7 flybys. Solid curve 
shows the intensity calculated using EUVAC model at SZA = 45°. Solid curve with symbols shows the limb intensity 
calculated using S2K model at SZA = 45°. Dashed curve shows the calculated intensity (using EUVAC model) at SZA = 
0°. Dash dotted curve shows the calculated intensity (using EUVAC model) at SZA = 0° and after reducing the e-C02 
cross sections forming CO(a'^n) by a factor of 2. Symbol represents the observed intensity of Cameron band and UV 
doublet band measured by Mariner 6 and 7 (Stewart, 1972). Observed values are shown for 2 orbits each of Mariner 6 
(for SZA = 27 and 0°; open and solid triangle, respectively) and Mariner 7 (for SZA = 44 and 0°; open and solid circle, 
respectively). 
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